Quantum walks: decoherence and coin-flipping games 
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We investigate the global chirality distribution of the quantum walk on the line when decoherence 
is introduced either through simultaneous measurements of the chirality and particle position, or as 
a result of broken links. The first mechanism drives the system towards a classical diffusive behavior. 
This is used to build new quantum games, similar to the spin-flip game. The second mechanism 
involves two different possibilities: (a) All the quantum walk links have the same probability of 
being broken, (b) Only the quantum walk links on a half-line are affected by random breakage. In 
case (a) the decoherence drives the system to a classical Markov process, whose master equation is 
equivalent to the dynamical equation of the quantum density matrix. This is not the case in (b) 
where the asymptotic global chirality distribution unexpectedly maintains some dependence with 
the initial condition. Explicit analytical equations are obtained for all cases. 

PACS numbers: 03.67.-a, 03.65.Ud, 02.50.Ga 



I. INTRODUCTION 

In the field of quantum computation, discoveries like 
the Shor and Grover algorithms [l|, [2[ have shown a supe- 
rior efficiency with regard to their classical equivalents. 
In this frame the quantum walk (QW) on the line Q, a 
natural generalization of the classical random walk, is re- 
ceiving permanent attention 4-6] . It has been used as the 
basis for optimal quantum search algorithms on sev- 
eral topologies and more recently it has been shown that 
any quantum algorithm can be reformulated in terms of 
a Q W algorithm 0, H[ , this means that the Q W is univer- 
sal in quantum computation. In particular the QW has 
the property to spread over the line quadratically faster 
than its classical analog. It remains a challenge to use 
this quantum property, as well as quantum parallelism 
and quantum entanglement, to increase the speed and 
efficiency of the algorithms. In this line of thought, it 
is important to probe further into the properties of the 
QW dynamics. 

On the other hand the spectacular development of 
techniques to manipulate atoms and photons has led to 
some experimental implementations of the QW as 
well as several implementation proposals [ll|. Here the 
obstacle of decoherence is present owing to imperfections 
and environmental noise. 

Several authors have studied the QW subjected to dif- 
ferent types of coin operators and sources of decoherence 
to analyze and verify the principles of quantum theory 
as well as the passage from the quantum to the classical 
world fl2l - fl6l |. These works have shown that the deco- 
herence transforms the quantum behavior into a classical 
diffusive behavior. In general, to obtain this conclusion 
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the authors focus their studies in the spreading of the 
wave-function. However decoherence may also drive the 
system to a more general behavior than the diffusive. 

m 

In recent work [18| the global chirality distribution 
(GCD) was defined as the distribution of the chirality in- 
dependently of position and the asymptotic behavior of 
the QW was investigated focusing on the GCD without 
decoherence. It was shown that the GCD has a station- 
ary long-time limit, a result usually linked to a Markovian 
process. 

In the present paper, we study the asymptotic behavior 
of the GCD when the system is subjected to two differ- 
ent sources of decoherence. In the first place, decoherence 
is introduced through simultaneous measurements of the 
chirality and particle position. This treatment allows us 
to connect the dynamics of the QW with certain aspects 
of quantum game theory [19j. In the second place, the 
decoherence is introduced through the influence of ran- 
domly broken links into the QW dynamics. This last 
mechanism may be relevant in experimental realizations 
of quantum computation based on Ising spin-1/2 chains 
in solid-state substrates [22j . 

In the next section we present the standard QW model 
and the GCD dynamics. In the third section we build the 
master equation for the GCD with periodic measurement. 
In the fourth section we present the coin-flipping games 
with the GCD. In the fifth section the GCD dynamics 
with broken links is studied, and in the last section we 
draw the conclusions. 



II. THE COINED QW AND THE GCD 

In this QW the walker moves (at discrete time steps 
t £ N) along a one-dimensional lattice of sites k £ Z, with 
a direction that depends on the state of the coin. The 
global Hilbert space of the system is the tensor product 
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H s (g) H c where H s is the Hilbert space associated to the 
motion on the line and H c is the chirality Hilbert space. 
The one-qubit "coin" subspace, H c , is spanned by two 
orthogonal states {\L), \R)}. The spatial subspace, H s , 
is spanned by the orthogonal set of position eigenstates, 
\k) . The evolution is generated by repeated application of 
a composite unitary operator U which implements a coin 
operation followed by a conditional shift in the position 
of the walker. This shift operation entangles the coin 
and position of the walker. In our case U depends on 
a parameter 9, with 9 € [0,7r/2] , that defines the bias 
of the coin toss (9 = ? for an unbiased, or Hadamard, 
coin). Then the unitary operator U{9) evolves the state 
in one time step as 



m+i)) = u{9)\*(t)). 



(i) 



The state of the total system at time t can be expressed 
as the spinor 
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a k (t) 
b k (t) 
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(2) 



where the upper (lower) component is associated to the 
left (right) chirality. Therefore, the QW is ruled by a 
unitary map, its standard form being [f| 



a k (t + 1) = a k+ i(t) cos 9 + b k+ i(t) sin( 
b k (t + 1) = Ok-i(i) sin 9 - bk-i(t) cos( 



(3) 



In references [23l . |24| it is shown how a unitary quan- 
tum mechanical evolution can be separated into Marko- 
vian and interference terms. This idea has been imple- 
mented recently in Ref. [18] for the global left and right 
chirality probabilities which are defined by 



P L (t)= £ |a fe (i)| 2 , 
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Pn(t)= E \b k {t)\\ 



(4) 
(5) 



k= — oo 



with P R (t) + P L (t) = 1. The pair formed by 



called the global chirality distribution (GCD). Starting 
from the original map Eq. © a quantum dynamical 
equation for the these distributions was obtained (ljj 



P L {t- 

P R (t- 



cos 2 9 sin 2 9 
sin 2 9 cos 2 9 



Re[Q(t)] sin26> 
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where 



Q(t)= £ a k (t)bl(t), 



(6) 



(7) 



with Re(z) indicating the real part of z. The two dimen- 
sional matrix in Eq. © can be interpreted as a transition 
probability matrix for a classical two dimensional random 
walk. On the other hand, it is clear that Q(t) accounts 
for the interference. When Q(t) vanishes the behavior 
of the GCD can be described as a classical Markovian 
process. 

It has been proven [H[ that Q(t), P L (t) and P R {t) 
have a long-time limit, and their values are determined 
by the initial conditions. Eq. © was solved in this limit 
defining 



U L = P L {t 
IIr = P R {t 

Qo = Q(t - 
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and substituting these values in Eq. ([5]). The stationary 
solution obtained for the GCD was 
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2Re(Q )/tan6» 
2Re(Q )/tan6» 



(9) 



Therefore, the dynamical evolution of the QW is unitary 
but the evolution of its GCD has an asymptotic value 
characteristic of a diffusive behavior. 

In this paper the initial state of the walker is assumed 
to be sharply localized at the line origin with arbitrary 
chirality 



l*(0)> = 



cos a 
expi/3 sin a 



|0>, 



(10) 



where a g [0, tt] and (3 € [0, 2ir) define a point on the 
unit three-dimensional Bloch sphere. Fixing the bias of 
the coin toss 9 = 7r/4, the analytical value of Qo was 
obtained in Ref. [25| following the method developed by 
Nayak and Vishwanath 26] 



cos2a + sin 2a (cos (3 + iV2 sin/3) 



Using this expression in Eq. 
limit of the GCD is 
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the asymptotic long-time 
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(12) 

If we assume that the parameters a and j3 verify tan 2a = 
— 1/ cos j3, the asymptotic solution is 11^ = 11^ = 1/2 and 
Re(Qo) = 0. In this case Eq. © approaches a Markov 
chain [27| with two states and the dynamics of the GDC 
turns into an example of dependent Bernoulli trials in 
which the probabilities of success or failure at each trial 
depend on the outcome of the previous trial. 

In Figs. [1] and [2] we present the probabilities 11^ 
and IIr, respectively, as functions of the parameters a 
and (3. In these figures we use 2000 time steps as an 
asymptotic temporal value. Zones with different tonali- 
ties of gray can be appreciated representing different oc- 
currence probabilities. The values of 11^ and Hr are ap- 
proximately restricted to the interval [0.3, 0.7] as a con- 
sequence of the dependence of Eq. (|12p on a and (3. 
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FIG. 1: The asymptotic probability II l as a function of di- 
mensionless angles using Eq. (|12[) . Seven zones of probability 
are given with different tonalities of gray. From the black to 
the white zone the intervals of probability are: (1) [0.29, 0.34], 
(2) [0.34,0.39], (3) [0.39,0.45], (4) [0.45,0.52], (5) [0.52,0.60], 
(6) [0.60,0.65], (7) [0.65,0.71]. In section IV this figure also 
represents Alice's strategic space for the game with one mea- 
surement 



FIG. 2: The asymptotic probability II_r as a function of di- 
mensionless angles using Eq. (|12l) . Seven zones of probability 
are given with different tonalities of gray. From the black to 
the white zone the intervals of probability are: (1) [0.29, 0.34], 
(2) [0.34,0.39], (3) [0.39,0.45], (4) [0.45,0.52], (5) [0.52,0.60], 
(6) [0.60,0.65], (7) [0.65,0.71]. In section IV this figure also 
represents Bob's strategic space for the game with one mea- 
surement. 



III. PERIODIC MEASUREMENTS IN THE 
CHIRALITY 



In this section we will consider a succession of unitary 
evolutions of the QW followed by a process of chirality- 
and-position measurement. Let us take the origin as the 
initial position of the random walker with the initial con- 
dition |W(0)) of Eq. (TO]), The position and chirality of 
the walker are jointly measured every T steps, with T 
sufficiently large in order to use the asymptotic results 
for the chirality Eq. JTJJ, i.e. T ~ 2000. The mea- 
surement of any observable of the system produces the 
collapse of the wave function into an eigenstate of the 
measured operator. Among the several alternatives for 
measuring chirality, we choose to measure it in such a 
way that the chirality is projected on the z direction by 
the Pauli operator tr z , whose eigenstates are the qubit 
states |L) = (J) and \R) = (°). The final position can be 
any position between k = —T and k — T but of course 
with different probabilities. Due to the fact that we are 
interested in the time-evolution of the chirality, we can 
rename the final position measured as the new position 
|0), this translation does not modify the chirality evolu- 
tion because this QW is homogeneous in position. Imme- 
diately after the measurement has been performed, the 
system is in state |L)|0) or |i?)|0) with probability LI^ or 



Hr respectively. 

The system evolves again unitarily during T time steps 
and then a second measurement is performed. Let us 
analyze this second evolution carefully. Supposing first 
that the system starts from the chirality state |L), after 
T time steps the new GCD is determined using Eq. (TT2"]) 
with a — 0, and the value obtained is (^). If instead 
we start from the chirality state \R) once again using 
Eq. (fT2"|) now with a = 7r/2, we obtain a GCD given 
by (*) . In the general case, we have a probability Hl of 
starting from the state \L) and II ^ of starting from the 
state therefore we must weigh the previous cases by 
their probability of occurrence. We obtain in this way 
the composite probability expression for the GCD after 
the second measurement 
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where 



p=l- 



2^2' 



q = 



2^2' 



(13) 

(14) 
(15) 



and we define P L (T) = U L and Pr(T) = II j 



Eq. (jT3|) is clearly a master equation and the global 
evolution of the system is a Markov process consisting of 
a succession of unitary evolutions followed by measure- 
ments. Such processes have the property that for any 
set of successive times (T, 2T, 3T, inT) the conditional 
probability at mT is uniquely determined by the value 
of stochastic variables at (m — 1)T and is not affected 
by any knowledge of the values at earlier times. In other 
words, the future state of the system depends only on 
the current state and not on the path of the process. In 
our case the measurement of the state of the system is a 
simple but extreme way of introducing decoherence that 
produces a loss of long-range memory. After performing 
m > 2 measurements the master equation becomes 



(16) 



and the probabilities as functions of the initial GCD evo- 
lution are 
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(17) 



where II l and 11^ are given by Eq. (|12p . After some sim- 
ple algebra, involving the calculation of the eigenvalues 
of the previous matrix in Eq. (fT6")) , the general solution 
of Eq. ([17]) is obtained 



where 



(18) 
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_ ( Pm Qrn \ 
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\ Qm Pm J 
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p m = \{l + {2 P -l) m - 1 } 
qm = \{l-{l-2qr- 1 } 



(19) 



In summary the global evolution of the system, 
throughout a time interval involving many measurement 
events, satisfies a master equation and therefore can 
be described as a Markovian process. The system has 
an unitary evolution only between consecutive measure- 
ments. For a sufficiently large number of measurements, 
Eq. ([15} implies that both Pl and P R , at first sight, 
tend to 1/2 independently of the initial conditions. In 
order to show this behavior graphically we present Fig. 
[3] where P R (2T) is calculated using the same color code 
used in Fig. [5J The values of P R (2T) are approximately 
restricted to the interval [0.39, 0.60]. This shows that the 
interval of probability [0.45, 0.52] has increased in com- 
parison to the case described by Figs. [Tj and [2j Addi- 
tionally, we also calculated the probabilities when three 
measurements are performed, that is P R (3T) with m = 3 
in Eq. (|T5|) . We verified that P R (3T) is approximately 
restricted to the interval [0.45,0.52], meaning that if we 
used once again the color code adopted before, the cor- 
responding figure would be an uniformly gray square not 
worthwhile presenting. 
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FIG. 3: The asymptotic probability Pr(2T) as a function of 
dimensionless angles using Eq. (|18[) with m — 2. Three zones 
of probability are given with different tonalities of gray. From 
the darker to the lighter zone the intervals of probability are: 
(1) [0.39,0.45], (2) [0.45,0.52], (3) [0.52,0.60]. In section IV 
this figure also represents Bob's strategic space for the game 
with two measurement. 



IV. COIN-FLIPPING GAMES 

Game theory has been used to explore the nature of 
quantum information. Initially quantum games were pro- 
posed as a quantum generalization of their classical coun- 
terparts but, due to the principles of quantum mechanics, 
new game possibilities have arisen within this scenario 
[19l - t21j . In the frame of the previous section, let us con- 
sider a simple quantum state flip game played between 
Alice and Bob. Initially the system is prepared in the 
arbitrary state |\&(0)) as in Eq. ([TO"]) where a is chosen 
by Alice and j3 is chosen by Bob, or vice versa. Next the 
QW dynamics evolves with Eq. (TTJ up to the time t > T. 
At that time the chirality of the system is measured. The 
players agree that the result of the measurement deter- 
mines who wins the game. For example, if the state is 
\L) Alice wins $1 (Bob loses $1) and if the state is \R) 
Bob wins $1 (Alice loses $1). This is a two-person, zero- 
sum game; the payoff given to Alice is the exact opposite 
of that awarded to Bob. If we suppose rational players, 
the election of their strategies is determined by the elec- 
tion of the initial condition for |W(0)). From Figs. [Tj 
and [H which show the dependence of probabilities II ^ 
and Tl R with a and /3, we can conclude that the first 
player to select one of the initial conditions has an ad- 
vantage. For example, if Alice chooses first, she needs 
to choose a £ (0.457T, 0.557r) (see Fig. [Tj because in this 
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way, independently of what Bob chooses, she makes sure 
that the probability to obtain \L) is 0.60 < ILr < 0.65 
in which case Bob almost certainly loses. On the con- 
trary, if Bob chooses first (see Fig. [5]) he needs to choose 
a E (0,0. In) or a E (0.9tt,tt) to obtain a probability 
0.60 < U R < 0.71 independently of /3. Therefore, the 
players must make their measurements obeying the rules 
of quantum mechanics and the initial conditions of the 
wave function determines their strategies. By looking at 
Figs. [1] and [H the players can plan their strategies and 
quickly conclude that this is not an equitable game; the 
first player has many winning strategies. 

A possible variant of the previous game is obtained 
when a second measurement at t > 2T is incorporated. 
In this new game, after the first measurement at t = T, 
the system has an unitary evolution until a second mea- 
surement is performed. Once again, Alice wins $1 and 
Bob loses $1 if the result of the last measurement is 
\L), independently of the result of the previous measure- 
ments. If the last measurement yields \R) Bob wins $1 
and Alice loses $1. Fig- El shows the probability Pr(2T) 
as a function of the initial conditions a and j3 for the 
new game. In this figure a large new area with proba- 
bility 0.45 < Pr(2T) < 0.52 appears. This area shows 
that this game is fairer because now the initial condi- 
tion /3 can be chosen in order to balance the Pl(2T) and 
Pr(2T) probabilities. If the number of measurements 
between unitary evolutions is increased, we obtain a suc- 
cession of more equitable games and in the infinite limit 
we have a fair game where Pl{-^ oo) = Pr{-^ oo) = 1/2 
for all initial conditions. 

Taking into account the inequity of the previous games, 
the players decide as a strategy to choose a and (3 at 
random, but obeying the constraint on 1 11/(0)) given by 
Eq. (TIT))) . With the aim to quantify the inequity of these 
games we now introduce the mean payoff of the game for 
each player. In the first game using Eq. (IT21 we can 
reinterpret II^ and Hr as the win density for Alice a a = 
TIl and Bob <tb = n# as functions of the parameters a 
and p. Then the probabilities of winning the game are 
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d/3 / a A {a,P) da = 1/2 



(20) 



d/3 / a B (a,/3) da =1 - ir A = 1/2, (21) 



for Alice and Bob respectively. The expected payoff for 
Alice is tta = %^a — $1ttb and the expected payoff for 
Bob is ttr = $l7Ts — $17ta; in this game tta = 7Ts = $0. 
In the other games, we define the win densities as in the 
first game but now using Eq. (TT5)) 
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1 - (1 - 2q) 
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2 ' 

2 ' 
(22) 



with m = 2,3, 4, .... 

The probabilities of winning the game are calculated 
using Eq. ((22]) in Eqs.([20l l2"T|) . Their values are n A = 
1/2, ttr = 1/2 again and the expected payoffs for Alice 
and Bob are tta = 7t_b = $0. Therefore, it is expected 
that they will tie on average after many games are played, 
provided they choose initial conditions at random. 



GCD DYNAMICS WITH BROKEN LINKS 
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FIG. 4: Situations that can arise at a site k of the line when 
there are: (a) no broken links (Eq. Q), (b) one link is broken 
fEqs. (|23l24[l h and (c) both continuous links are broken (Eq. 
(I25p ). The upper white (lower black) circles indicate the left 
(right) chirality amplitude. The arrows indicate the direction 
of the probability flux in one time step. At each circle two 
arrows point in and two point out in order to conserve the 
total probability flux. 



Let us now consider the breaking of links as a mech- 
anism to introduce decoherence in the QW [28]. This 
model was presented for the first time in Ref [13j. Sup- 
pose that, at time t, a given site k may have one or both 
of the links connecting it to its neighboring sites bro- 
ken. If site k has no broken links, as shown in Fig. 0] 
a, then the evolution law Eq. ([3]) is applied. When one 
or both links at site k are opened there can be no trans- 
lation across the broken link and the evolution must be 
modified accordingly. If the link to the left of site k is 
broken, as shown in Fig. 0] b, the upper component of 
the chirality at k receives probability flux from k + 1. In 
order to preserve the flux, the outgoing probability flux 
from the upper component at k must be diverted to the 
lower component at the same site. The corresponding 



6 



transformation on the chirality components is therefore 

a k (t + 1) = a k+1 (t) cos 9 + b k+1 (t) s'm9, 
b k (t+l) = a k (t) sin6» +b k (t) cos 9. (23) 

If the broken link is to the right of k, the situation is 
similar and the transformation is 

&k(t + 1) = a k(t) cos 9 — b k (t) sin 9, 

b k (t + 1) = a k -i(t) sm6 - b k -i(t) cos9. (24) 

Finally, if site k is isolated, as in Fig. 0] c, the unitary 
operation is followed by a chirality exchange, 



a k (t 
b k (t 



1) = a k {t) cos ( 
1) = a k (t) sin t 



■ b k (t) sin( 
b k {t) cost 



(25) 



Hence, there are four situations that can arise at a site 
k: (a) no link is broken, (b) the link on the left side is 
broken, (c) the link on the right side is broken, and (d) 
both links are broken. The global evolution of the sys- 
tem depends on the application of one of the four maps 
in each spatial position, where each map has associated 
its correlative operator. Thus the dynamical evolution of 
the GCD can be characterized by four unitary operators 
which will be called in generic terms U n with n = 0, 1, 2, 3 
where Uq will be related to the operator U of Eq. ([1}. 
A statistical description can be obtained combining the 
operators into a single evolution equation with the ap- 
propriate statistical weights. 

To implement the algorithm of the quantum walk with 
broken links we proceed as follows. At each time step t, 
the state of the links in the line is defined. Each link has a 
probability r of breaking in a given time step, being r the 
only parameter in the model. Therefore the probability 
that a given site has no adjacent broken link is 



r = (l-r) 2 , 
that it has a left or right broken link is 

r\ = r 2 — r(l — r), 
and that it is isolated is 

r 3 = r 2 . 



(26) 



(27) 



(28) 



These probabilities are appropriate statistical weights 
since they satisfy 



To + n + r 2 + r 3 = 1. 



(29) 



In order to study the chirality dynamics of the QW we 
introduce the reduced density operator 



Pc = tr/?, 



(30) 



where the operator 'tr' is the partial trace taken over the 
positions and p is density matrix of the quantum system 



(31) 



Using the wave function given by Eq. ([2]) and its normal- 
ization properties, the reduced density matrix obtained 
in Ref. [29| is 



Pc 



( P L (t) Q(t) 

\Q(ty p R (t) 



(32) 



This matrix depends explicitly on the GCD and its time 
evolution is given by the time evolution of the GCD. The 
dynamical equation for pJt) is obtained using the density 
matrix formalism I15l 1 1 71. 1301 



Pc (t + l) = Y,A nPc (t)Al 



(33) 



71=0 



where a set of Kraus operators ({^4n} with n = 0, 1, 2, 3) 
has been introduced in order to simulate the decoher- 
ence produced by the broken-links process. The Kraus 
operators are defined by 



A n — \ 



(34) 



and, to preserve the trace of the density matrix, they 
satisfy 



N 



A-n A^ — I, 



(35) 



where / is the identity matrix. 

We want to study the quantum evolution of Eq. (|33p , 
but this task is not trivial since we do not know the set 
of operators {A n } explicitly. However, we can infer some 
knowledge about the quantum evolution of the GCD if 
we follow the procedure developed in section II to ob- 
tain Eqs. ((6|) from Eq. (|3]). Starting from the maps 
Eqs. (|3l23l24l25j) the corresponding dynamical equations 
for the GCD are built. Next these equations are weighted 
with the probabilities r n to finally obtain the master 
equation for the averaged GCD 



P L (t+l) 
P R {t + l) 



cos 2 9 sin 2 9 
sin 2 9 cos 2 9 



PUt) 
P R (t) 



(36) 



This equation has no dependence with the parameter r 
and, in opposition to Eq. ([6]), has a negligible interfer- 
ence term Q(t) due to the random phases introduced by 
the lack of coherence between the terms of the sum that 
define the interference in Eq. ([7]). Therefore, if Q(t) is 
negligible it is clear that Eq. ( |3"3"|) is equivalent to Eq. 
(1361) . and we could use this equivalence to obtain infor- 
mation about the operators {A„}. 

The stationary solution of Eq. (|36l) is = = 1/2 
independently of the initial conditions and of 9. We im- 
plemented a numerical code with the four maps Eqs. 
PI23I24I25[) . weighted with the values r , ri, r 2 , r 3 re- 
spectively, and have verified the asymptotic value I4l — 
1/2 for several initial conditions and several values of 9 
and r. Thus, the theoretical result given by Eq. (|36[) has 
been numerically confirmed. 



7 



0.78 



0.76 



I1 L 0.74 




0.72 



0.78 



0.76 



I1 L 0.74 




0.72 



0.72 0.74 

Re(Q )+1/2 



0.78 



FIG. 5: The asymptotic probability II .l as a function of the 
probability r with 8 = 7r/4 and the initial condition 775 (J) • 
The stars are the values obtained from the numerical calcu- 
lation and the thin line is a polynomial adjustment. The 
calculation has been made using the map Eq. ((3| for the left 
side of the evolution and the maps Eqs. (|3I23I24I25|) weighted 
with the values given by Eqs. (|26I27I28[1 for the right side of 
the evolution. 



FIG. 6: The asymptotic probability FIl as a function of the 
dimensionless term Re(Qo) + 1/2, with 9 = n/4 (see Eq. ((9J) 
and the initial condition Q). The stars are the values 
obtained from the numerical calculation and the thin line is 
a linear adjustment with slope 1. The calculation has been 
made using the map Eq. (O for the left side of the evolution 
and the maps Eqs. (|3I23I24I25[) weighted with the values given 
by Eqs. (|26l27l28p for the right side of the evolution. 



The usual QW dynamics, Eq. ([3]), has the property of 
spreading over the line quadratically faster than its clas- 
sical analog and after a long time its wave- function is very 
dispersed on the line. This situation makes it interesting 
to study the decoherence when only some segment on the 
line is affected by the breakage of links. Such decoherence 
mechanisms may be relevant as a model of experimental 
realizations in quantum computation [ToL ITT1 . [22j as well 
as in any problem related to the transfer of quantum in- 
formation through space using QW [1, Q. In order to 
investigate this type of decoherence, we propose a com- 
posed non-homogeneous QW with two different dynamics 
starting form the origin: for the left half-line the evolu- 
tion is determined by Eq. and for the right half-line it 
is determined by the set of equations Eqs. (|3I23I24I25[) . 
weighted with the values ro, r\, r%, respectively. The 
stochastic behavior of the GCD , for an ensemble of these 
QWs, will depend, for the left side, on the usual QW and, 
for the right side, on the QW with broken links. There- 
fore the right half-line chirality distribution satisfies Eq. 
(fMf and the left half-line chirality distribution satisfies 
Eq. ©. From this reasoning, it is clear that the GCD 
is the sum of the above right and left partial distribu- 
tions and it satisfies Eq. (jfi]), where the interference term 
emerges from the amplitudes of the wave function belong- 
ing to the left side of the evolution. However, it should 



be noted that the incorporation of the interference term 
of the right side of the evolution does not modify this 
equation because it remains negligible. 

The evolution of the system ceases to be unitary owing 
to decoherence, however in spite of the stochastic nature 
of the system Re(Qo) 7^ and Eq. © maintains the de- 
pendence between the asymptotic values Hl and Re(Qo)- 
The system loses its integrability and the explicit value of 
Ul is to be calculated numerically. We have performed 
this calculation, as a function of the probability r with 
9 = 7r/4 using the original maps Eqs. (|3I23I24I25[) with 
an ensemble of 100 dynamical trajectories and 2000 time 
steps. Figure [5] shows that Hl has a parameter depen- 
dence on r. Figure [6] shows the linear dependence be- 
tween Ul and Re(Qo) + 1/2, with slope 1 for 9 = n/4, 
verifying the stationary solution given by Eq. ©. There- 
fore, we conclude that there is a perfect agreement be- 
tween the theoretical stochastic approach presented in 
this section and the numerical calculation using the orig- 
inal maps. We also found from the numerical calculation 
that the asymptotic distribution 11^ depends on the ini- 
tial conditions. Figure [7] presents the numerical calcula- 
tion of the distribution II l , with 3000 time steps and the 
initial condition given by Eq. (fTDj) . 

Comparing Fig. [7] with Fig. [TJ it is possible to distin- 
guish several similarities, such as the number of light and 
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FIG. 7: The asymptotic probability II i, as a function of 
dimensionless angles. The calculation has been made us- 
ing the map Eq. <[3j) for the left side of the evolution and 
the maps Eqs. ()3I23I24I25|) weighted with the values given 
by Eqs. (|26I27I28|) for the right side of the evolution, with 
r = 0.3. Seven zones of probability are given with different 
tonalities of gray, from the black to the white zone the inter- 
vals of probability are: (1) [0.738,0.746], (2) [0.746,0.755], 
(3) [0.755,0.763], (4) [0.763,0.771], (5) [0.771,0.780], (6) 
[0.780,0.789], (7) [0.789,0.797]. 

dark spots in both figures. These similarities, show that 
the asymptotic behavior depends on the initial condi- 
tion. In other words, this type of decoherence erases the 
unitary behavior of the system but does not completely 
erase its dependence on the initial condition. However, 
the range of values of the probability is now restricted to 
~ [0.7, 0.8] while in the Fig. Q]this range was ~ [0.3, 0.7]. 
Such a decrease in the range shows that the trend of the 
system is to go to a single asymptotic distribution that 
will only depend on the parameter r, as seen in Fig. [5] 

VI. CONCLUSIONS 

Decoherence in quantum systems like the QW has 
been extensively studied. Analytical and numerical re- 
sults on the effect of different kinds of noise have shown 
that quantum properties are highly sensitive to random 
events. In this paper, we studied the asymptotic behav- 
ior of the QW on the line when it is subjected to different 
sources of decoherence. We focused on the dynamics of 



the GCD when decoherence is introduced in two different 
ways. In the first case decoherence is introduced through 
periodic measurements of position and chirality. It has 
been shown that the evolution of the GCD, in a time 
scale involving many measurements, can be described as 
a Markovian process and its master equation has been 
obtained analytically. This master equation allows us to 
build new quantum games, similar to the usual spin-flip 
game, where as a novelty the players perform measure- 
ments on the QW system. These games are characterized 
by a continuous space of strategies, and the selection of 
an initial condition determines the particular strategy 
chosen by the player. These games may be used as a 
tool to study quantum algorithms subjected to external 
decoherence, as in the extreme case of measurement. 

In the second case, the decoherence is introduced 
through the influence of randomly broken links affect- 
ing the QW dynamics. Two different variants of this 
second case have been investigated. These are: (a) All 
the QW links have a certain probability of being broken, 
the GCD is described by a Markovian process. Its mas- 
ter equation is obtained analytically and it is shown to 
be equivalent to the dynamical equation of the density 
matrix. It was shown that the dynamics is independent 
both of the breakage probability and of the coin bias of 
the QW. (b) Only the links on one of the half-lines (say 
the right one) are affected by the possibility of break- 
age. In this case the dynamical evolution of the GCD 
is not described by a Markovian process but it has an 
asymptotic value that depends on the initial conditions 
and the breakage probability. Here the interference term 
Qo is not negligible but an asymptotic stationary relation 
between Qo and the GCD has been analytically found. 
The behavior of this composite QW is at first sight un- 
expected since usually decoherence destroys the unitary 
correlation, providing a route towards a classical-like be- 
havior described by a Markov process. To understand 
this behavior suppose a localized initial condition far to 
the left of the origin, situated thus in a zone where the 
evolution is determined by the usual QW map. As time 
passes, the wave function spreads in both directions and, 
due to the dynamical evolution of the map, the extreme 
left branch of the wave function will never be influenced 
by the existence of decoherence introduced into the right 
branch of the wave function. In this way, the wave func- 
tion keeps some information about the initial condition. 
We conclude that the effect of decoherence of the sec- 
ond kind studied in this work does not necessarily trans- 
form the quantum system into a dissipative one such as 
a Markov process. In more general terms, the mere pres- 
ence of noise does not assure the transition of a quantum 
system to classicality. 

We acknowledge the comments made by V. Micen- 
macher and the support from PEDECIBA and ANIL 
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